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Abstract 



Liouville field theory on a sphere is considered. We explicitly derive a differential equation 
for four-point correlation functions with one degenerate field V_rnb. We introduce and study 
■ also a class of four-point conformal blocks which can be calculated exactly and represented by 

. finite dimensional integrals of elliptic theta-functions for arbitrary intermediate dimension. We 

study also the bootstrap equations for these conformal blocks and derive integral representations 
for corresponding four-point correlation functions. A relation between the one-point correlation 
^ . function of a primary field on a torus and a special four-point correlation function on a sphere is 

^ [ proposed. 

1 Introduction 

Liouville field theory has attracted a lot of attention since the discovery of its application in the 
quantization of the string in non-critical dimension in 1981 [1]. It gives an important example of a 
non-rational unitary conformal field theory (CFT) with central charge c > 25. From another side 
any correlation function in minimal models i. e. in rational CFT's with central charge c < 1 can be 
obtained from the appropriate correlation functions in LFT by analytic continuatioiu. 

Liouville field theory is defined on a two-dimensional surface with metric Qab by the local Lagrangian 
density 

C = ^g-'^da^dk^ + ixe''"^ + ^R^, (1.1) 



^At this point some care is needed (see Ref [2]). 
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where R is associated curvature. This theory is conformally invariant if the couphng constant h is 
related with the background charge Q as 

Q = b + \. (1.2) 
The symmetry algebra of this conformal field theory is the Virasoro algebra 

[L^, L„] = (m - n)Lm+n + ^f*^^^ ~ -m (1-3) 

with central charge 

CL = l + &Q^. (1.4) 

Primary fields in this theory, which are associated with exponential fields e^"''', have conformal 
dimensions 

A(a) = a(Q-a). (1.5) 

The fields and Vq^a have the same conformal dimensions and represent the same primary field, i.e. 
they are proportional to each other: 

V^ = R{a)VQ.^, (1.6) 

with the function 

(7r/i7(6^))(Q-^")/^ 7(2&a-&^) 
^(«) - ^2 - 2a/h + 1/62) ' 

known as the reflection amplitude. Here and later we use the notation 

7(x) = r(x)/r(i -x). (1.7) 

Since Liouville field theory is conformal invariant it can be solved exactly at least in principle. In 
particular it means that one can find explicitly correlation functions of all local fields. Following the 
ideas of the operator product expansion (OPE) it is sufficient to find the structure constants of the 
OPE in order to find all correlation functions in the theory. In the case of Liouville field theory this 
problem simplifies drastically due to conformal invariance of the theory [3]. Namely, the operator 
product expansion of two primary fields V^^ and has the form 

-1 poo Q 2 

V^,{z)V^M = - / Vq^,p(0) + . . . dR (1.8) 

Here by . . . we denote the contribution of the so-called descendant fields i.e. fields which can be 
obtained by the action of the negative part of the Virasoro algebra (11.31) on the primary fields; this 
contribution is universal and completely determined by conformal invariance [3]. The part of the 
operator product expansion which is not fixed by conformal invariance is the set of structure constants 
of primary fields C^f^j • They can be found from the condition of associativity of the operator algebra 
[4]. It is better to write down the three-point correlation function 

(Vai (21)^.2(^2)^3 (2:3)) = |^^^|2Ai+2A2-2A3|2_^3|2Ai+2A3-2A2|223P^2+2A3-2Ai' (^-9) 

which is trivially related to the structure constant C"l^^ as 

Caia2 =C(ai,a2,Q-a3). (1.10) 
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The constant C{ai,a2,a^) introduced by (11. 9p was proposed in [5,6] and passed a lot of crucial tests 
in [7] (more rigorous proof can be found in Ref. [8]) 



C(ai, 02,^3)= 7r/i7(6^)6 



[,2-252 



T(6)T(2«i)T(2a2)T(2a3) 



where a 
functional relations 



(1.11) 
satisfying 



;i.i2) 



T{a — Q)T{a — 2ai)T(a — 2a2)T{a — 2^3) ' 
ai + 0:2 + «3 and T{x) is an entire function, self dual with respect to b —>■ b 

T{x + b) =^{bx)b'-^''''T{x), 

T{x + r^) = 7(rV)6'''"'"-^T(x), 

which for general real values of the parameter b^ have a unique solution given by the integral flA.ip 
with normalization T{Q/2) = 1. The analytic part of the operator product expansion (11.81) i. e. the 
contribution of descendant fields contains more non-triviality. Despite the fact that this contribution 
is determined completely by conformal invariance its determination represents a tedious problem. This 
problem was attacked in Ref.s [9-11] where analytical properties of the so-called four-point conformal 
block were studied. The four-point conformal block is defined as a contribution of the particular 
primary field together with its descendants in four-point correlation function. Namely, due to (II. 8p 
the four-point correlation function can be represented as a sum over intermediate states 



{Va^{Zi, Zi)Va2iz2, Z2)Vas{z3, Z3)Va^{Zi, Z^)) = JJ 



7.1 27ii • 



i<j 



x^j c(a^,a2,^ + ^p)c{^-^P,as, 



012 ^3 
«! 0:4 



X 



dP, (1.13) 



here jij are standard combinations of conformal dimensions [3] and x is the anharmonic ratio for the 
four points Zj: 

[l.U) 



X 



^13^24 



The contour of integration C in (I1.13P goes along the real axis with possible deformation due to the 



presence of discrete terms (see Ref. [7] for details). The conformal block 



OL2 03 



X ) sums up all the 



intermediate descendant states of a given primary one with conformal dimension A 



+ P2 in the 



operator product expansion (11.131) . T-channel conformal block usually is represented graphically as 



(^2,0) 



5i 



Ol2 "3 



X 




P2 + ^ 




(as, 00) 



;i.l5) 



("4, 1) 



Unfortunately this function is not known in a closed form, but it can be found as a power series 
expansion. The most efficient way to do that was proposed in Ref. [11]. According to [11] we introduce 
instead of usual coordinate x a new coordinate r 
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which is the coordinate on the upper half plane (Im(r) > 0). Here K{x) is the elliptic integral of the 
first kind 



Kix] 



1 



dt 



The elliptic conformal block ( 



2 Jo v^t(l-t)(l -xt)' 
q ) is defined as follows 



;i.i7) 



Ci2 a-i 



X ) = (16g) a; 



q;2 as 



:i.i8) 



where g = e™'^ and ^3(0') is the theta constant (see definition in appendix^])- The function q 

satisfies a nice recursive relation (see Ref. [11] and also Ref. [12]) which leads, in particular, to an 
effective algorithm for calculation of its expansion in power series of q 



a2 0:3 



\ai 0:4/ 



1.19) 



Due to the much better convergence of this sum than in the usual variable x the representation of 
the conformal block in terms of the elliptic variable q is more convenient for numerical studies. The 

coefficients Sj^p^ ( ) are rational functions of 

■'f \^ai on J 

^(i)r«2 «3^ _ 32&^((ai - - {a, - ^)\{a, - f )^ - [a, - 



APW + (62 + 1)2 



;i.20) 



In principle, using the recursive relations any coefficient Sjp alj expansion fll.lQp can be found 

explicitly, but they are rather cumbersome for large values of L (even for L = 2 the answer is very 
complicated). In order to find the four-point correlation function ( 11.13^ one has to perform the 
integration over the intermediate momentum P in f ll.lSp . Altogether it represents a complicated 
numerical problem, which was studied for particular cases in [7, 12,13]. In some cases the four-point 
correlation function can be found in explicit form. In particular, if one of the fields is degenerate (see 
below) the correlation function satisfies a differential equation of Fuchsian type with solutions having 
explicit integral representations. 

This paper consists of two parts. The first part (section [2]) is devoted to the derivation of the explicit 
form of the differential equation for the four-point correlation functions in Liouville field theory. In 
the second part of this paper (section [3]) we discover a family of conformal blocks which have explicit 
integral representations for arbitrary value of intermediate momentum P. We study the bootstrap 
conditions for these conformal blocks and derive the four-point correlation functions. We propose also 
a relation between one-point correlation function on a torus and special four-point correlation function 
on a sphere in LFT. In appendixes we collect some formulae used in this paper and give proofs of some 
statements. 

The content of this paper is closely related with the subject where the seminal contribution of 
Alyosha Zamolodchikov has a great importance. The authors had a chance to know him personally 
and could estimate his scientific virtuosity, originality and creativity. We dedicate this paper to the 
memory of Alyosha. 
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2 Differential equation for four-point correlation function 

Among the primary fields Va the so-called degenerate fields are of special interest. They form a closed 
subalgebra of the operator algebra and are characterized by the parameters 

mb n 

a = amn = — Y ~ 2b' ^'^''^^ 

Any degenerate field V^„„ has a null- vector in its Verma module at level (m + l){n+ 1). In particular 
the primary field V_b has a null- vector at the second level 

{L\ + b'^L_2)V_i=0. (2.2) 

As a consequence of the null-vector condition (12.21) the correlation function which contains this degen- 
erate field and an arbitrary number of general primary fields! 

{V_.iz)V^M)---VaA^n)) (2.3) 

satisfies the second order partial differential equation [3] 



{V.iz)V^M)---VaA^n))=0. (2.4) 



When the number of external points Zk equals 3 this equation can be reduced to the ordinary Riemann 
differential equation (which is equivalent to the hypergeometric equation) 



, k=l 



A(Qfc) 
{Z - Zkf 




{V_i,{z)V^,{z{)V^,{z2)V^,{z;)) = 0, 

(2.5) 



where 



A ( ) and A12 = A(ai) + A(a2) - A(a3) etc. 



It means that the four-point correlation function with one degenerate field V_h can be expressed 

in terms of the hypergeometric function^. A more general four-point correlation function with one 
degenerate field V"_mf> 

(y_^{z,z)Vo,MVa,{l)Vo,A^)) (2.6) 



satisfies an ordinary differential equations of order m + 1 in both variables z and ^ as a consequence 
of the null- vector condition. A procedure to find differential operator which "kills" the correlation 
function (12.61) from the null- vector condition is straitforward and was described in Ref. [3]. However 
this procedure becomes tedious for large values of the parameter m and it is difficult to construct such 
a differential operator explicitly using this method. Instead in Ref.s [14, 15] an explicit representation 



^Here and later we will drop sometimes the dependence on the antiholomorphic variables z. 

^One has to remember that a differential equation similar to (|2.5p but with opposite chirality z ^ z is also valid and 
the four-point correlation function is a bilinear combination MijHi{z)Hj[z) where Hi[z) and Hj{z) are solutions of the 
corresponding differential equations. 
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for the correlation function (12.61) was obtained. Namely, the correlation function (12.61) can be expressed 
in terms of a 2m- dimensional Coulomb integral 



J^(A,S,C|x) 




n 

fc=i 



^ \2A\j. 1 |2B|. „|2C 



rf^ti . . . 



(2.7) 



with parameters 

A = 6 (a - 2«i - g + m6/2) , B = h{a-2a2-Q + mb/2) , C = 6 (Q + m6/2 - a) 
and the normalization constant Qm{ai, 02,0^3) is 



fi„(ai, 02,03) = (-Ti-yu)™ vr/i7(6 )6' 



I,2\r,2-262 



(Q-ci-m6/2) 



T'(-m6) nLiTr(2«^) 



T(a - g - f ) riLi '^(c, - 2afc + f ) 



(2.7a) 



(2.7b) 



here a = 01 + 02 + 03. Integration in (12. 7p over each two-dimensional variable goes over the 
plane. This integral is convergent for some values of the parameters o^ and b, otherwise it should 
understood by means of analytic continuation. The result (12.71) was obtained from arguments which 
did not involve the explicit construction of the differential equation. However, for several important 
purposes one needs the differential operator for the four-point correlation function in explicit form. 

In this paper we find the differential equation for the correlation function (12.61) . Partially we follow 
the logic of paper [16]. Let us consider the five-point correlation function with one degenerate field 



{V^iz)V^MVa,im,ioo)V^A^)). 



{21 



This correlation function satisfies a second order partial differential equation. For future purposes it 
is convenient to define the function '^{u\q) as 



{V_^{z)V^,{Q)V^,{l)V^,{^)V^,{x)) = z^{z 



, iz{z-l)iz-x)Y QAuf ' , , , , , 



2A(a4) 1 
(X(X - 1)) 3 +12 



0'i(o) 



where the variable u is related with the variables z and x as 

z — x 

TT j a:(z-l) dt 



u 



(2.10) 



and 61 (u) is the Jacobi theta function (see definition in appendix El). We note that (12.101) maps a 
two-sheeted covering of the sphere with four marked points x, 0, 1 and 00 onto a torus T with periods 
TT and 7rr where r is given by (11.161) . Marked points x, 0, 1 and 00 are mapped to the points 0, |, ^ 
and T^^^ on a torus (see fig. [1]). The function ^{u\q) defined by eq (12.91) satisfies the non-stationary 
Schrodinger equation with doubly periodic potential 



dl 



where the potential U(u) is given by 



4? 

V{u) + —d^ 

TTO^ 



U(m) = V(u) + 6Y 



(2.11) 



(2.12) 
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f 

Figure 1: Transformation (12.101) maps a two-sheeted covering of the sphere with four marked points 
onto a torus which is represented by the parallelogram with the periods vr and vrr. Points x, 0, 1 and 
oo are mapped to the points 0, |, ^ and ^'^y^ on a torus. 



with 

4 

Y{u) = J2sjisj + '^)pi^-^j) (2-13) 
i=i 

and parameters Sk are related with parameters as 

". = |-^(.^. + i). (2.14) 
In (12.131) p{u) is the WeierstraB elliptic function with periods vr and vrr defined by the infinite sum 

^ ^ ' (2.15) 



{u — Txk — -ktIY {irk + tttIY J ' 



and u;,- are the half periods 



71 TTT 

ui = -, UJ2 = — , tus = u;i + u;2, = 0. 



(2.16) 



We choose the w-independent term 6Y in (12.121) in such a way that the expansion of the potential 
U(M|a;) at the origin has the fornfl 



l]{u\x) 



S4(S4 + 1) ^^Vr('=+^)(r 



k=l 



.2k 



-U 



(2.17) 



where 
with 



TT 



2A: 



^^^^^ ^ (-l)'=+^(x-l)wf^P,,(x)+xwf )Pfe(l-x)-(-l)'=-V(x-l)wf x'-2Pfc(l/x) 

(2.17a) 



w 



(fc) 



+ i) + £ii£i±l) 



(2.17b) 



*This u-independent term is SY = — X]fc=i ^kisk + l)p{u>k) 
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In eq f l2.17al) Pk (x) are polynomials of degree k — 2 





— X 




— X 


^4(1 


— X 


^5(1 


— X 



4 
3 



^(2 + 13a; + 2x2), (2.18) 
5 

64 

-—(1 + 30a; + 30x^ + 0;='), 
35 



These polynomials can be obtained as an expansion of the square of the elliptic sine (see appendix | 
for connection between the WeierstraB function and the elliptic sine function). Namely 



sn2(t|v/^) = f;^±iii^P. (2.19) 



k=l 



k\ 



The polynomials Pfc(l — x) defined by eq fl2.18p can be rapidly calculated using the recursive formula 



m(2m - l)P.„+i(l - 3^) = 3x^ _ ^^,7^'^ _ ^^,2 ^fc(l ' x)P^+i,fc(l - x) (2.20) 



k 

with initial condition 



Pl{l-X) = -Ul + ^]- P2{l-X) 



It is important also to have in mind the expansion 



^n-\t\v^) = l + \{x+i) + Y. ^'^a (2-21) 

k=l 

where 

Gk{x) = -±-(^{-lf-\x - l)Pfc(x) + xPfc(l - x) - (-l)^-ix(x - l)x'=-2Pfc(l/x)). 

Let us try to find a solution to the differential equation (12. lip in the fornj^ 

^{u\q) = u''+' (^(r) + ^_i(r)M2 + ^_2{r)u^ + ■■■), (2.22) 
then (12.111) will be equivalent to the semi-infinite matrix differential equation [16] 

fc=i / 



^As we will see later the function ^'(t) defines the four-point correlation function of the fields Va,. with k = 1,2,3 and 
the field V^^_j^. The functions '^^^{t) define the four-point correlation function of the fields Va^ and some descendants 
of the field V'" i . 
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J+ = 







1 





^(r) = 




) 




-2s4 - 5 
















1 




) 









-S4-| 


^) 












^) 






/ 



(2.24) 



Or explicitly 

(S4 + ^)vi>_i + 4^9.VI/ = 0, 

2 T[b^ 
(2s4 + 5)^_2 + + W^^\r)^< = 0, 



We see that if the parameter S4 in eq fl2.14p takes the valued 

3 



S4 = —m — 



(2.26) 



then the infinite chain of equations (12.231) has a finite sub-chain which can be written in the form 
(12.231) . but with finite (m + 1) x (m + 1) matrices J+ and jj^: 



J_ 



/ 


. 








o\ 


m 


. 














. 


. . 2(m-l) 










. 





m 


0/ 



/ 





1 




























1 

0/ 



(2.27) 



From this finite chain of equations we can conclude that the function \E'(r) in eq_^ (l2.22l) satisfies a 
differential equation of the order (m + 1). This equation can be written in the forn]^ 



d: 



{tri+l 



+ 



m(m + l)(m + 2) 
6 



im— 1 



mir) = 0. 



(2.28) 



All functions W^''^ (r) will enter in the differential equation (I2.28P with positive integer coefficients. It 
is reasonable to change variable r — > x in fl2.28p using projective invar iance. Namely, the differential 
equation (12.280 is invariant under the change of variables 



r 



{duj\ 



1 



W^^\u;) + -{u;,t}, W^''\t) 



/duo\^ 
' \d^) 



iy(*=)( 



(2.29) 



00 



for k > 2, 



®The general term below the diagonal in J_ in (|2.24p is equal to — ri(s4 + n + i). 

''it corresponds to the situation Oi4 = ^ — ^ and hence in the operator product expansion V_j_ (2)1^04 (x) appears 

the degenerate field V_ . 

2 

^More precisely, in this case we meet the situation when one of the solutions to the differential equation (|2.1ip 
contains logarithmic terms in the variable u (the so-called resonance case). The condition that these logarithmic terms 
are cancelled is equivalent to the condition that the finite-dimensional subchain of differential equations (j2.23p with 
matrices J+ and J_ given by eqs (|2.27p is satisfied. 

(2 \ k 
2L|- W^W(r). 
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where {f,z} is the Schwartz derivative, which is defined as {f,z} = f"'/f' — ?)/2{f" / fY . Choosing 

K{l-x) dr 



K(x) 



dx 4:x{x — l)K'^{x) 



and using that in this case 



{x,t} 



+ 



2x2(1 -a;) 2{x-iyx 2x{l - x) 



(2.30) 



(2.31) 



we obtain that the function \l/(x) satisfies the same differential equation which can be easily derived 
from the system 02.231) with matrices J_|_ and J_ given by fl2.27p and the substitution ^ 



m+l 



m{m + l)(m + 2) 
6 



W^^\x)d^-' + ...]<il{x) = 0, 



(2.32) 



but with currents W^^\x) given b}0 

wS'^Pfc(x) 



x] 



+ 



wf^Pfe(l 



\Nl^^X^-^Pk{l/x) 



x'^ll — x)'^ ^ {x — l)^x'^' ^ 



1 — X 



ifc-1 



(2.32a) 



where 



w 



(2) 



W 



(k) 



i2k-2 



\ ^ 2 Ah 



L2k 



with 



4^= - 1; n = 3 ■ 4'^-^ 



(2.32b) 



if > 2 



L (2.32c) 

and polynomials Pfc(l — x) defined by eqs (12.181) . fl2.19p . We note that the function ^'(a;) in fl2.32p is 
related with the four-point correlation function with one degenerate field V_rnb{x) as (we write down 
only the chiral part) 



x-'^'\x - l)-^'\V_^{x)V^Am.Al)V..{oo)) 



(2.33) 



We have checked the first seven differential operators obtained by the above method with those obtained 
by the "brute-force method" and found a complete agreement. We list the first few examples of the 
differential equations in the appendix [Bl We see that the coefficients before all terms in the differential 
operator f l2.32p are positive integer numbers which can be easily obtained from the matrix equation 
f l2.23p with matices J+ and J_ given by f l2.27p and have a combinatorial structure. 

It is instructive to notice that the coefficients W'^'^^x) in the differential operator (12. 320 transform 

as 

W^^\x) {-lfW^^\x), (2.34) 

under the substitution b ^ ib, ^ 1 — A^, which corresponds to the differential operator killing 
the four-point correlation function with one degenerate field in the complementary theory with central 
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Here we are going back and use the parameters Afc — ak{Q — ak) instead of the parameters Sk defined by (|2.14p . 



10 



charge c = 26 — and with additional fields with conformal dimensions = 1 — A^. Transformation 
f l2.34l) refiects a very simple property of the differential operator of type fl2.32p with arbitrary functions 
(currents) W^''\x) [17]. Namely, if ipi, . . . y'ipm+i are some (m + 1) linearly independent solutions to 
(12:321) then Xi, . . . , Xm+i; where Xk is given by the Wronskian 



Xk 



i'l ■■■ i^k-l 'ipk+l ■ ■ ■ V'm+l 

dtpi . . . di)k-i di)k+i . . . di^rn+i 



(2.35) 



are solutions to the differential equation fl2.32p with currents W^''\x) — >■ {—l)''W^''\x). It means 
that the corresponding correlation functions in both theories are simply related. Another interesting 
property of the differential operator fl2.32p under the transformation (12.341) is that it transforms up to 
the sign (—1)™-+^ to the hermitian conjugate operator. 

Let us mention also that we obtain an important example of a differential operator with a unitarized 
monodromy matrix. In general, if we consider the differential equation (12.320 with currents given by 
(I2.32ap but now with arbitrary numbers w^*^^ and polynomials Pk{x) of degree (fc — 2), then it gives the 
most general Fuchsian differential equation with three singular points 0, 1 and oo with singularities 
defined by Without loss of generality we can assume that the polynomials Pk{x) are normalized 

as 

Pfc(x) = 1 + . . . (2.36) 

The subleading coefficients in polynomials (I2.36P define so called accessory parameters which do not 
affect on singular behavior but of course contribute to the monodromy. The problem is to tune them in 
such a way that monodromy is unitarized. With a requirement of that type we deal with the solution 
of the si[n\ Toda on a sphere with three singular points [18]. The unitary monodromy allows one to 
build a bilinear combination of holomorphic and antiholomorphic solutions which is a single-valued 
function on a sphere with three punctures. In paper [18] it was proposed that the solution to this 

(k) 

problem is unique in some special domain of the parameters . It means that all coefficients in 
(I2.36P i.e. accessory parameters are fixed. However this statement is not yet proved (at least to our 
knowledge) nevertheless it looks reasonable from a physical point of view. The fixing of the accessory 
parameters is a very transcendental problem and its solution is far from being done, but in some 
cases one can guess an answer. Here we have found a class of differential equations with unitarized 
monodromy which is parameterized by four numbers Ai, A2, A3 and b and the accessory parameters 
or what is the same the polynomials Pk{x) given by (12.180 . It follows trivially from the fact that in 
this case the solution to the corresponding differential equation represents the holomorphic part of the 
four-point correlation function (12.70 which is obviously a single-valued function because it is given by 
a multiple integral over the plane. This class of differential operators can be a starting point for the 
numerical study of the accessory parameters problem. 

All results of this section can be trivially rewritten for the correlation function with degenerate 
field V-EL, but for general degenerate fields V^„„ with m 7^ and n 7^ this construction is more 
complicated and we will not discuss it here (see also section H]) . 



"'^^By gauge transformation ^'(a;) x"{x— l)^'i'{x) one can always set the coefficient before the term with subleading 
derivative to zero. 
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3 Integrable potentials and conformal blocks 

In the previous section we constructed the differential operator for the four-point correlation function 
(12. 6p which has an explicit integral representation. In this section we consider the differential equation 
(12.111) and for special choice of the parameters find an explicit integral representation for its solution. 
By the transformation 

^(M|g) -^x^^(x- l)^2^(M|g) (3.1) 

where 

Ai = ^ (2 (ai - Q/2f - ia2 - Q/2f - (a^ - Q/2f) , 

A2 = ^ (2 (a2 - Q/2f - (ai - Q/2f - {a, - Q/2f) , 
the differential equation (12.111) can be transformed to the generalized Lame heat equation 

Ai 

{-dl + V(w)) ^{u\q) = —drm{u\q), (3.2) 

where the potential V(m) is given by (12.130 . For integer Sj the potential (12.130 is known to be a 
finite-gap potential and called the Treibich-Verdier potential [19]. The differential equation (13.21) was 
studied in particular cases in Ref.s [20-22] (see also Ref. [23] where interesting relation with Baxter's 
Q— operator for the eight- vertex model was discovered). We propose, that if all parameters Sk in eqs 
(|27[3D - (I27[D take the values 

Sk = mk + -^, (3.3) 

then the differential equation (13.21) is integrable in the sense that the general solution to (13. 2p can 
obtained from the general solution to the heat equation: 

dl^{u\q) + -%d^^{u\q) = (3.4) 
no 

by an appropriate integral transformation. For example let Si = S2 = S3 = and S4 = 1, then the 
general solution to the equation 

4? 

{-dl + 2piu)) vl/(«|g) = :^d^^{u\q), (3.5) 
can be obtained from the general solution \I^o(w|g) of the heat equation (13. 4p as follows: 

Jo \Q'iiO)3 J E[u)E[v) 
where we introduced the notation 

9i(m) 

E{u) = (3.7) 

The proof of (13.60 and of more general relations can be found in the appendix O In principle, the 
integration in (13. 6p can go between any two zeroes of Qi{v), for example between and vr or between 
and vrr. It is instructive to consider the limit h ^ 00 and take in (13. 6p 

^o(« + h\\q) = g^e-^("+''"), (3.8) 
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with A being fixed at 6 — oo. In this case the integral in (13.61) is governed by the saddle point z/, which 
is solution to the equation 

and \I'('u|g) has a limit (up to irrelevant factors): 

m^) - ^7^^-^". (3.10) 

which is the solution to the stationary Lame equation with energy p(z/). So, the solution (13. 6p can be 
viewed as a "quantization" of the solution to the stationary equation (13.101) . More suprising is that 
equation (13.21) is integrable for values of the parameters Sk which vanish in the limit 6 — ^ oo. Let us 
consider for example the case Si = S2 = S3 = and S4 = ^, then the general solution to 

2 / 2 \ \ Ai 
< + P + ^(")J ^(^1^) = (3-11) 

is given by the integral transform 

.H,.e;,o)i-^)/- (1^)* *,..2H,)*„ ,3., 

where \E'o('u|g) is the solution to (13.41) . In the general case (13. 3p solution is likely to be given by an 
integral of dimension 

= 5f + rii + ^2 + ^3 + ?T,4, (3.13) 
where g is the number of gaps for the classical potential [24] 

g = ]^ (^2maxmfc, 1 + m - (1 + (-l)"") (minm^ + \^ , (3.14) 

here m = "^mk- We justified this hypothesis in many cases where an explicit solution has been 
constructed. We believe that this conjecture is true, however for general integer values of and 
an explicit integral representation for the solution at present is not known. In many particular cases 
four-point conformal blocks (see below) which can be viewed as a limiting value of the solution to (13. 2p 
are given by the integrals of elliptic functions. Some of them are listed in appendix [El 

Although we assume that equation (13.21) is integrable for all values of the parameters Sk given by 
(13. 3p the solution can be very non-trivial. All explicit examples known to us show that the solution 
depends irregularly on the parameters Sk and it is difficult to obtain closed expression for all values 
(13. 3p . Below we consider the case Si = S2 = S3 = and S4 = m, which corresponds to the five point 
correlation function 



{V_±{z, z)Vi (2^-i)<, (x, x)K,(0)K,(l)K,(oo)) 
z{z-l) 



x{x — 1) 



"""^l ill f .M4(m+lf- ''"-'lf^'' |ei(^|g)|^ 

\Z — \X{X — 1)\ ^ ^' ^ 2 2^m[U,U), (3.15) 

|0;(O|g)|^+3 
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where 

^2 4' 

The function ^rn(w, u) satisfies the differential equatior''^^' 



Ai 



(-dl + m{m + l)p{u)) *rri(M, u) = u) 

whose general solution is given by the m-dimensional integral (see appendix [C]) 



(3.16) 



(3.17) 



^{u\q) 








n 

fc=i 



em 



n 

Kj 



<^l(v^ 



0'i(O) 



-b2 



n ll^V^^X ^o{u + h\\q)dv,...dv^, (3.18) 
E{u)E{vk) 



where v = Vi + ■ ■ ■ + Vm and \E'o('w|q') is again some solution to the heat equation 03.41) . Apparently, the 
space of solutions to (13.171) is infinite dimensional and spanned by the integrals (13.181) with \E'o(u|g) 
taken to be equal to 

v|/±(„|g) = qP'e^^'''P\ (3.19) 

Each solution \l/p(u|g) corresponds to a particular five-point conformal block. In this paper we consider 
four-point conformal blocks which can be obtained from the five-point conformal blocks in the limit 
M — »• 0. More precisely, there are two types of four-point conformal blocks which can be obtained 
from the solution (13.181) in this limit. This is related with the fact that in the operator product of 
the degenerate field VLj_ with any other fields in the correlation function (13.151) there appears only 
two primary fields togetlier with their descendants. In particular, in the operator product of the field 
V_± with the field Vi (2m-i)b we obtain two terms (we used the refiection relation a ^ Q — a in the 

26 2b 4 

second term) 



V_J_Vl (2r, 
2!- 2b 



-1)6 



V 



(2m-l)6 
4 



+ 



(2m+3)b 



(3.20) 



The first term in the r.h.s. of (I3.20p corresponds to the main asymptotic of the solution (I3.18P at 
M which is \l/(M|g) ~ u~"^ and gives the conformal block 



def ^ / 2 



_ b 

4 

(2m-l)b 



(3.21) 



where S)p{. . . \q) is defined by (11.181) . Fortunately, in order to obtain the conformal block (13.211) one 
has to take instead of \l/o in (13.181) one of the simplest solutions (13.191) (it does not matter which one, 
because both of hem have the same asymptotic at m — ^ 0). Having in mind definition (13.151) of the 
function "^rn{u,u), we obtain the integral representation for the conformal block (I3.2ip 



nl^\q) = K 




^2bP{ui+-+Um) 



llE{ukr''l[\E{ 



Uj — M,- 



dui . . . dUr, 







with normalization constant Nm given b\[^ 



N„ 



n 



fc=i 



k=i 



r 1 



2 



r 1 + 



(2m+l-fc)fe^ 



m 

(m-H)&2 



(1 - 1) nr=i r (1 + ^ + ihP) r (1 + ^ - ihP) 



(3.22) 



(3.22a) 



^^It satisfies also the antiholomorphic differential equation similar to (|3.17|) but with u u and t —> t. 
^^Normalization constant (|3.22ap and more general normalization constants were calculated in the appendix [Pj 
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We have performed an expansion of the integral fl3.22p in series at g — and compared it with a known 
expansion following from Alyosha Zamolodchikov's recursion formula for the conformal block [9] and 
find complete agreement up to high orders. 

It is interesting that not only the conformal block has a simple expression for this special choice of 
the external conformal dimensions. The product of structure constants which enters in the definition 
of the four-point correlation function fll.lSp simplifies drastically and is equal tc0 



4 '2 4' 2 \2 '2 4' 2 4 



iTjj,'y{b )h 



2\ j2-2f)2 



2 2b2 

X 



4T(6)2T (11^) T(|) ™ / / ,..2S 

X _^A_J_A_ll^ - ^) -< {-O'P - ^) ■ (3-23) 

As a consequence the integral over the intermediate momentum P of the product of two structure 
constants f l3.23p with the modulus squared of the conformal block (13. 220 which in fact gives the four- 
point correlation function can be performed analytically. We note that when the parameters Sk of the 
fields take values fl2.14p the dependence of the product of structure constants on the momentum P 
can be expressed as a product of 2A^ 7-functions where N is given by (13.130 depending linearly on 
this parameter. This fact supports our conjecture that in this case the solution is given by an integral 
of dimension N. The integration over the momentum P is rather non-trivial mainly because the 
contour of integration is deformed as shown on figure O This deformation of the contour is prescribed 
by the condition that the four-point correlation function is single-valued. Surprisingly, the result of 
integration over the momentum P is given by a multiple integral over the torus T with periods vr and 

7TT 

\qr'dnXP\r)dm{-P\T*) 



'c nr=i sin {i^iihP + ^)) sin [ixiibP - ^)) 

= A^(^lm(r)) /.../ \\£{uk,Ukr^' W£{ui - Uj.Ui - Uj)-''" d?ui . . .d\m, (3.24) 



T T 



k=l i<j 



where 



7r 7r 



k=l i<j 







and 



S{u, u) = E{u)E{u)e~^^ (3.24b) 

_ (-l)"^m! 2^'"-^sin'"(^) 

nfc=ism(^) sm( ^ 2 ) 

A phenomenon of this type is well known in the minimal models of CFT where the sum of squared 
modulus of conformal blocks which given by the contour Coulomb integrals taken with appropriate 



^''in order to simplify this product and obtain l|3.23p it is convenient to use the double argument formula for the 
T-function (|X2|l . 
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' f + lb-' 




Figure 2: Contour of integration C in (13.241) . This picture is drawn with the assumption that 6^ <^ 1. 
Otherwise the contour has to be deformed. 



coefficients can be always expressed in terms of integrals over the plane and this property trivially 
ensure single- valuedness of this correlation function [25]. Here we meet exactly the same situation 
but with the integral over intermediate the momentum P. This property was verified numerically 
with high accuracy (results of numerical calculations can be found on-line at [26]). We intend to give 
analytical proof in separate publication. The function S{u,u) given by (I3.24bl) is nothing but the 
exponent of the Green function of the Laplace-Beltrami operator on a torus 



u 



Trimr 



(3.25) 



Finally we obtain the expression for the four-point correlation function 



{V (2^-1)6 (x,x)v;,(o)K(i)K,(oo)) = Cr 





e2{q)BA{q) 


4m+2m(m+l)b^ 








\x{x - 


- 1)|5 + - 





X 




T T 



k=l 



Uk,Uk 



\mb' 



Y\_^{Ui — Uj, Ui — Uj) ^ Sui 



.d\^, (3.26) 



where 



C 



(^) 



7r"'m! 



7r/i7(6 )b 



2 2b2 



n 



7(-¥ + 



The right hand side of (I3.26P up to a trivial factor looks like a Coulomb gas representation of the 
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one-point correlation function of the operator V^mb' in LFT with parameter 6' = on a toru£j. We 
propose that this equahty holds for general fields (i.e. for arbitrary not necessary integer values of m). 
Let us define function T{a,b\q) related to the one-point correlation function 



{Va)r = TT(q 



,2Lo-^^2Lo-^ 



(3.27) 



in Liouville field theory with cosmological constant /i and coupling constant 6 on a torus with modular 
parameter r as 



2\i2-2b'^ 



\ri{r)\ 



-4A(q) 



(3.28) 



where rfij-) = {i202{q)dz{q)0i{q))^3 is Dedekind eta-function. We define also the function S{a, b\q) which 
is related to the four-point correlation function in LFT on sphere as 



1,2-262 



- + — - 



' \x{x - l)p^(") {V^{x,x)V^{Q)V^{l)V^{^)). (3.29) 



The correspondence between the one-point toric and the four-point spheric correlation functions states 
that 

f- b 
\a - -)V2,— 

V2. 



V2 



(3.30) 



Sia,b\q)=^(^{ 
where 'R{a,b) is given by 

We propose to give a proof of this relation for arbitrary values of the parameter a in a different 
publication. We note here that the relation between the conformal dimension of the operator on a 
torus (we denote it as A^) and the conformal dimension of the operator on a sphere (we denote it as 
As) in fl3.30l) has the form 

1 . 1 36^ ^ ^ 
As = -Ar + - + — . 3.31 

2 4 16 ^ ' 

In conclusion we note that all results of this section can be considered in the theory with c < 1 
(6^ < 0). We did not discuss this problem here. We stress only that our results can be useful for 
the calculation of four-point amplitudes in Liouville gravity (see for example [27,28]). Not going 
into details, let us consider the analytical continuation of Liouville four-point correlation function 
(13.261) to the region 6^ < 0, where it can be considered as a correlation function in matter theory with 
central charge 26 — c^. The corresponding correlation function of dressing fields (fields with dimensions 
Afc = 1 — Afc, where A^ is conformal dimension of matter field) in Liouville theory will be 



(\/(2^(x,x)V;,(0)K,(l)K,(oo)), 



(3.32) 



i. e. exactly the correlation function which is complementary to (13.261) in a sense that it is described 
by the second term in the operator product (13.201) . This second term defines the conformal block 



dof ^ / 2 

= rip 



b Q 
4 2 

(2m+3)b Q 

4 2 



_ b 
4 

b 

4 



(3.33) 



-'^^On a torus the one-point correlation function of the operator V-mV can be screened by m screening fields, because 
curvature term does not contribute to the total charge in this case. 
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which corresponds to the solution of (13.171) with asymptotic \1/ ~ u"^'^^. This solution however has 
more complicated expression in terms of elliptic functions. For example the conformal block Ti,[^\q) 
has the form 

7i(^)(g) =iVfi r e^''''^E{ufp\u)du, (3.34) 



where the normalization constant A^^i is given by 



(52_i)(62_2)r(f -^6P)r(f + ^6P)- 



The same is true in the more general situation when all Sk is eq (13. 2p are integer (s^ = m^). Both 
correlation functions in matter and Liouville theory in this case have integral representations comple- 
mentary to each other. We plan to study the corresponding four-point amplitudes in Liouville gravity 
elsewhere. 



4 Concluding remarks 

The differential equation defined in section [2] has solution given by the integral (12.71) . Similar integrals 
appear in different applications in CFT, in supersymmetric CFT [29] in CFT's having higher spin 
symmetry like Toda field theory [18,30], in different perturbed models of CFT etc. The results of the 
section [2] give us a correspondence between the four parametric family of differential equations and 
the integrals (12.71) . These differential equations permit effectively to calculate the integrals (12. 7p for 
all values of x. 

Here we considered the family of the operators VI mb or V-il. For general case of the operator 
V_rnb_ri_ thc differential equation is of order (m + l)(n + 1) and can be written in the form fl2.32p but 

2 2b 

with functions W^^\x) which are not anymore linear functions of the conformal dimensions A(afc) of 
the fields V^^. We note however, that the solution to this equation can be represented as a bilinear 
superposition of the solutions to differential equations of orders (m + 1) and (n + 1). In particular in 
the simplest case m = n = 1 the solution to the fourth order differential equation can be represented 
as a bilinear combination of hypergeometric functions. The correlation function for the field V_ 
= V_q{x, x) with three arbitrary fields up to a constant can be represented by the Integra.^ 



6 _ J_ 
2 26 



{V_q{x,x)V^, (0) V^, (1) V^, (cx))) ^ Ixl'"^'^ \x - l|'"^«/(pi,p2,P3,a;) (4.1) 

2 

where 

I{PUP2,P3,X) = j \t - 1|'^^ |t - \U - 1|'^^ 1^ _ ^\^P'z \t _ ^2^ ^2^ (42) 

here pi = b{a — 2«i — Q/2), p2 = b{a — 2«2 — <5/2), Ps = b{3Q/2 — a); with a = ai + 02 = «3 and 
P'i/Pi = b~^- 

It is convenient to use the notations pu = P1+P2, Pus = P1+P2+P3, then a useful integral relation 
for I{pi,P2,P3) takes the form: 

nPi,P2,P3,x) = C\xf+'^'' \x-lf-^'^''Hip,,p2,P3,x) (4.3) 



"'^^This integral has a singularity at t —> u and should be understood as a bilinear combination of contour integrals 
(see [31]). 
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where the constant 

7(-Pl23) 



C = -(l + b' 



i{-pih{-P2h{-P3) 

and 

H{p,,p2,P3,x) =J |tr'^-' \t - ir'^'-' \t - x\-'''~' \uf^'^-' \u - \u - x\^^'^-' \t-u\^ dHdPu. 

(4.4) 

The last integral can be easily rewritten in terms of contour integrals which can be expressed in terms 
of hypergeometric functions. Four conformal blocks in this correlation function can be classified by 
the parameter a[ of the intermediate primary field, which takes values: a'^ = ai — Q/2, 012 = ai + Q/2, 
ctg = ai — Q/2 + 6, 0:4 = 0:4 — Q/2 + 1/6. We normalize our conformal blocks G^'. by the condition: 

It is convenient to introduce the functions: 

J^i{pi,P2,P3,x) = F(l +p3,2 + pi23,2 +pi3,a;), 
J^2iPuP2,P3,x) = x~^~P^''F{l+p2,-pu-pn,x), 

where F(a, b, c, z) is hypergeometric function, and to define 

G^,{x) = a;i+Pi3+"iQ(i _ (4.6) 

then: 

Qa'Ax) = J='2{pi,P2,P3,x)J^l{-p[-l, -P2, -P3,a;)-y^ ^^^^f^ (Pl " 1, P2 , PS, S^) J^l (-Pl , -^2, -p's, x) 

(4.7) 

The function Qa'.^{x) can be derived from Qa{{.x) by the substitution pi —* —p[ 

Gadx) = ^^^P^^^^^^ ^^^^ :Fi{pi,P2,P3,x)J^i{-p[ - 1, -P2, -P'3,X)- 
' {P2 + P2) 

(1 -K23) a+Pislj,^^^ _ i,p,,^3,^)j-^(_p;, (4.8) 



{P2 + P2) 

and 

P3(Pl23+K23) n / \ T/ 1 A-r-^' ' 'A 
7Z , / = ^2 Pi - l,P2,P3,a;P^2 -Pl,-P2,-P3.a; - 

Pl3(l -Pl3)(l +P13) 

- ^2{Pi,P2,P3,x)T2{-p'i - 1, -P2, -P3,a;). (4.9) 

We note that in case when all four fields are V-q/2 the correlation function (14.11) was studied in [32] 
where it was expressed in terms of the hypergeometric functions. 

In this paper we did not consider the application of elliptic conformal blocks and correlation 
functions to physical problems and two-dimensional quantum gravity. We note only that for 6^ < 
the correlation function (I3.26P has a logarithmic behavior at coinciding points and probably can be used 
in physically interesting logarithmic CFT's. Correlation function (12. 9p appears also in studying SU{2) 
WZNW models [33] as well as in their non-compact versions associated with s/(2) and [34]. In 
these cases the variable z plays role of an isotopic variable. The results of this paper can be applied to 
studying of correlation functions in these models. We propose to consider the application of correlation 
functions derived in this paper to the two-dimensional quantum gravity in a future pubhcation. 



19 



Acknowledgment 



This work was supported, in part, by RBRF-CNRS grant PICS-09-02-91064. The work of A. L. was 
supported by DOE grant DE-FG02-96ER40949, by RBRF grant 07-02-00799-a, by the Russian Min- 
istry of Science and Technology under the Scientific Schools grant 3472.2008.2 and by the RAS program 
"Elementary particles and the fundamental nuclear physics". A.L. and E.O. thank the Laboratoire de 
Physique Theorique et Astroparticules Universite Montpellier II for hospitality. 



Appendix A. Useful formulae 

T(a;) function. This function is defined by integral representation: 

2 „;„T 2 (Q 



log Tfa;! 



~ dt 
T 



Q Y , sinh^ - x) I 

a; e -r^ 

2 / sinh y sinh ^ 



Double argument formula for the T-function: 

24x(x-S)+i 



T(2a;) 



T(|)T(^) 

Shift formula for the T-function (here we put a lower index 6 to T): 

b\ _ T46)T,(|) 



b \ V2 



Theta function. The theta function QAuW) defined b}o 



(A.l) 



(A.2) 



(A.3) 



^(n+|)2g(2n+l)m 



(A.4) 



is solution to the differential equation 



Ai 



a^ei(M|r) drQiiu\T)=0. 

IT 

It satisfies the following quasi-periodicity relations 

ei(M + 7r) = -Qi{u), 

ei(M + 7rr) = -g-^e-2^"ei(M), 

and transforms as follows under the action of the modular group 



ei(M|r + 1) = i20i(M|r), 



(ir) 2e 6i('u|r). 



(A.5) 



(A.6) 



(A.7) 



^^During this paper we drop sometimes the dependence of the function 8i(w|r) on r. 



20 



For some purposes it is useful to define the function E{u) 

E{u) = = sin(M) + 4sin3(M)g2 + 12 sin=^(M)g^ + O(g^). (A.8) 
The following differential relation will be useful in the appendix O 

3A - ^-3.) 1^ = 0. (A.9) 

TT / E[u)E[v) 

This relation can be also rewritten as 
E"{u + v) ^ E"{u) ^ E"{v) ^ E{u + v) fE'ju) ^ E'{v)\ ^ ^E^E^ _&m ^ ^^.10) 



E{u + v) E{u) E{v) E{u + v) \E{u) E{v) J E{u) E{v) 0;(O) 
Other theta-functions can be expressed through 9i(m) as 

e2H = ei(M + |), e3(n) = e™g3ei(n + | + ^), e4(w) = -ze™gi0i(n + ^). (A.ll) 

We define also theta constants 

^2(g) = 02(0), eM = ©3(0), e,{q) = 64(0), (A.12) 

and functions 

E,(u)^E(u), Mu)^^, E:M^^, E,,)^^^. (A,13) 

Weierstrafi function and Lame equation. The WeierstraB function is defined by the infinite sum 
(12.151) and can be expressed through the second logarithmic derivative of the theta- function 9i(m) as 

It is related to the elliptic sine function as 



The WeierstraB function has the expansion at the origin 

1 , 5-2 ^ 2 , fi'3 4 



= - + ^u' + + Oiu'), (A.16) 
20 28 



where the numbers g2 and also known as invariants are given by 

4 f2K{xV ^ 



3^ = ^y-^j (a; -^+1), 

4 f2Kix)V , 
^73 = 7^ —] {x + l){x-2){2x-l). 
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27 V TT 



(A.17) 



We give also the double argument formula for the WeierstraB function 



k=l 



where ujk are half periods given by (12. 161) . The Lame equation 

{-dl + 2p{u))^{u) = -p{v)^{u) 

has a solution 

, E(u + v) ^2i^u 
^ ' E{u)E{v) 
which trivially leads to the following identity 

E(u + v) 



E(u)E(v) 



e'liv) ie;"(o)\ e{u + v 



Qiiv) 3 e;(o) J E{u)E{vy 



(A.18) 

(A.19) 
(A.20) 

(A.21) 



Appendix B. Covariant differential operators 

Here we give the first few examples of covariant differential operators defined in section [2) 



^(n+i) ^ _^ + l)(n + 2) ^^^(2)(^^^)g« 



-1 



The coefficient before the derivative of the order (n — A;) are some graded differential polynomials of 
order (fc + 1) in the functions (currents) W'^'^\x) . . . W^''^^\xj^ Some coefficients of these polynomials 
are fixed from the condition that currents transform in appropriate way under the action of the 
diffeomorphism group 

W(^\x)-.(^^yW^'\u;) + ^{u;,x}, ^ (^) V^H for k > 2. ^^'^^ 

Namely, up to order 5 all coefficients are fixed from this condition (up to total normalization) 

D(2) =52^H^(2)(a;)^ (B.3) 

= d^ + AW^^\x)d + 2dW^^\x) + W^^\x), (B.4) 
= d^ + 10W^^\x)d^ + {mdW^^\x) + QW^'^\x))d+ 



2)(5) = 20W^^\x)d'^ + {30dW^^\x) + 21W^'^\x))d^+ 

+ {6AW^^\xy + 18d'^W^^\x) + 21dW^^\x) + 8W^^\x))d+ 
+ {32d{W^^\x)Y + Ad^W^^\x) +6d^W^''\x) +AdW^''\x) +A8W^^\z)W^''\x) + W^''\x)). (B.6) 



^^Currents W^''\x) have weight k and derivative dx has degree 1. 
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First ambiguity appears at level 6, because one cannot distinguish between the functions W^^\x) and 
(W^^^x))"^ using only their transformation properties. Both of them transforms like a tensor, so in 
order to separate them some other symmetries should be taken into account. Our analysis in section 
[2] gives the differential operator S)^^-* and all other operators. For example S)^^-* equals 



+ i259W^^\xf + 63d^W^^\x) + SAdW^'^^x) + 36W^^\x)jd^+ 
+ (259d{W^^^ + 28d^W^^^ (x) + ASd'^W'-^^ (x) + 36dW^*^ (z) + AAOW^^^ {x)W'-^^ (x) + lOW^^^ (x)) d+ 
+ (5d^W^'^^ (x) + lOd^W^^^ (x) + lOd'^W^^^ (x) + 5dW'-^^ (x) + 155W^'^^ {x)d'^W^'^^ (x) + 130 {dW^^^ {x)Y + 
+ 220d{W^^\x)W'-^\x)) + 225{W^^\x)f + 100W^'^\x)W^^\x) + m{W^^\x))^ + W^'^\x)). (B.7) 



As we see from eqs (IB.3P - flB.7P the form of these differential operator is canonical in a sense that their 



coefficients are integer numbers which can be easily obtained from the WZW equation fl2.23p with 
finite matrices J+ and J_ given by fl2.27p . 

Appendix C. Integral representation for solutions to the gen- 
eralized Lame heat equation 

We consider differential equation 

4 



4? 

^{u\q) = —dr^{u\q). {C.l] 



-dl + ^ Sk{Sk + l)p{u - Uk) 
k=l 

As was argued in the section [3] this equation is integrable for values of the parameters = ruk + ^ 
We consider for simplicity the case Si = S2 = S3 = and S4 = 1. We have 



4? 

[-dl + 2piu)] ^iu\q) = :^d^^iu\q). (C.2) 



Let us try to find a solution in a form 



^{u\q) = I P{^\<D^^^fo{u + h\\q)dv, (C.3) 



where we choose the integration limits in such a way that surface terms can be neglectecuj. Applying 
the differential operator (—9^ + 2p{u) — -^dr) to the integral flC.3P and using identities flA.QD and 



f lA.2ip we find that the function p{v\q) should satisfy two equations 

1 d,p{v\q) _ @'^{v) &i{v) Ai d,p{v\q) _ 1 <$>"M 

62 p{v\q) ei(t;)' Q^{v) Tx\? p{v\q) 3 0;(O) 



(C.4) 
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It will be justified below, that it can be done. 
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with solution 



p{v\q) 



e;(o)^ 



(,2 



(C.5) 



So the contour of integration goes between any two zeroes of Qi{v). 

Generahzation to the case Si = S2 = S3 = and S4 = m is very straitforward. We have 



[-dl + m{m + l)p{u)] 

The solution has the form 



l62 



^{u\q) 




n 

k=l 



euo)3 



E{u + Vk) 
E{u)E{vk) J 



n 



Ai 



&l{Vi - Vj) 



(C.6) 



euo)3 



'^o{u + b'^v\q)dvi...dVm, (C.7) 



where v = vi + ■ ■ ■ + Vm- It can be verified again with the help of relations (\A.9\i and (]A.21|) . Acting 
with the differential operator [—d^ + m{m + l)p{u) — ■:^dr] on (lC.7p . using the trivial identity 



\Eiu)Eivk) J 



and integrating by part (we assume that integration limits in (\C7h allow to do that) we obtain an 
integrand which is proportional to 



( E"ivk) l er(o) 
-iKEivk) 3e;(o) 



E'{vk) f E'{u + Vk) E'{u) 



k=l 



E{vk) \E{u + Vk) E{u) 



2E 

Kj 



E'{u + Vi) E'{u)\ (E\u + Vj) E\u 
E{u + Vi) ~ E{u] 



Kj 



E"{v, - vA 



E{u + Vj 

1 0;"(o) 



E{vi-Vj) 3e;(o) 



E{u) 
+ m(m — 1) 



-25: 

E{u) 



E'{u + Vk) E'{vk-Vj) 
E{u + Vk) E{vk - Vj) 

' E"{u) i e;"(o) ' 
E{v) ^ 3 e;(o) 



(C.9) 



Expression ( jC.91) can rewritten as a sum of two terms 



A 



[m 



^^lE{u + Vk) E{u) 



E"{vu) _ + ( E'{u) E'{vu) \ 

E{vk) E{u + Vk) \E{u) E{vk)J 



+ 



+ 2 



E'{u)E'ivk) er(o) 
e;(o)J 



(C.IO) 



and 



!E"(v, 



i<j 



Vj) E"{u + Vi) E'\u + Vj 



L — fj) E{u + Wj) + V 



- 2 



^'(^* - Vj) [E'{u + Wi) ^'(m + Vj] 



E{vi — Vj) \ E{u + Vi) E{u + Vj) 
_ ^ E\u + v,)E\u + v,) _ &[iOy 

E{u + Vi) E{u + Vj) e;(o). 



(C.ll) 
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In both sums (IC.lOp and (IC.lip each term is equal to zero identically due to (lA.lOp . 

One can easily find an integral representation for the solution to (IC.ip in the more general case 
Si = S2 = S3 = and S4 = m + ||^. We give it without a proof (which is however straitforward and 
based on using identities f[01) . flA.lOp and flA.2ip ) 



^{u\q) = (0;(O)) 



2n 1-12-, 




n 



\ ve;(o)l 



E{U + ffc) 



n 



E{u)E{vi 



n 



©l(fi - Vj) 



e;(o). 



n 

Kj 



0'i(O) 



n I "^l!!. I ^o(m + b\ + 2t;'|g) rf"^; dV, (C.12) 



0'i(O) 



where = t;i + ■ ■ ■ + Vm, v' = v[ + ■ ■ ■ + v'^, dJ^v = dvi . . . dvm and d^v' = dv[ . . . dv'^. The integrand in 
flC.12p has a singularity at Vi —>■ Vj and should be correctly defined as a principal value similar to what 
was done in Ref. [35]. One can obtain a solution to fIC.ip with arbitrary Sk = rn + ^ (not necessary 
S4) and other sj = just by the substitution E{u + v^) Ek{u + Vk), E{u + f^) Ek{u + v[) and 
E{u) —>■ Ek{u) in flC.12p . where Ek{u) are defined by flA.lSp with identification 



Si-^Ei{u), Si^E2{u), S2^E4{u) 



S3 E^iu) 



(C.13) 



Appendix D. Normalization integrals 

Here we give some integrals which can be used as normalization factors for conformal blocks. We 
introduce the notations: 



-^"(m) = JJ^ |2sin(uj - , d'^u = dui...dun, g = -iP' /2. (D.l) 

The first normalization integral that appears in the normalization of conformal blocks has the form: 

In{a,p,g) = ne2^"'= {2sin{uk)f D^iuf^d^u (D.2) 
7 — 1 







fc=i 



This integral can be calculated exactly and is equal to: 

in^a^P.g) \yV{l + g)V{l + l-g + jg + ip)V{l+-^-g + jg~ip) 
For the our purposes we need the values of this integral for a = nb"^, g = —6^/2 and p = bP 



f.\ r (1 - f ) r(i + f + i6P)r(i + f - ibP) 

The modification of this integral that also appears in the applications is: 



(D.3) 



(D.4) 



In{nb\ bP, -b^) = Lf f[ e'''-f'"'= (2 sm{uk))"'' cot (m,) Dniu)-'' d^'u. (D.5) 
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This integral is equal to I„ = J„(n6^, bP, — 6^/2)/i„(P), where 



(2n + l-j)|- 



Another normalization integral appears, when we consider conformal blocks given by integrals from 
—TIT to TTT avoiding zero clockwise for pure imaginary r and take the limit r — > ioo (we did not 
consider them in this paper, but they still are important): 

oo oo ^ 

Ua,p,g)= !...! J]e2'P"'^-(2coshK))'^D„(iM)23rf"M (D.7) 

-oo -oo '^--^ 

This integral can be reduced to the Selberg integral and is equal to: 

II \ r(i + jg) r(-f - g + jg + ip)r(-f -g + w- W ,„„, 
J..(a,P,9) = n^(^ 2r(-. + g-,g) ■ <°-«) 

In particular 

• • }ir(l-f) 2r(-(2n + l-j)f) * ' 

The integral analogous to (ID. 50 in this case is 

oo oo ^ 

J„(n62, 6P, -672) = j-j tl e''^^"' (2 cosh(Mfc))"''' tanh (uk) Dniiuy^' (Tu (D.IO) 



k=l 



which is equal to: Jn{a,p,g) = i"/in (-P) Jn{'nb'^,bP, —6^/2). 

For some conformal blocks considered in the appendix [E] (eq flE.SP for k = 2) normalization factors 
are defined by the integral: 

MP) = y!.^e^^^"/^(2sin(«,))(^""^)/'' (2cos(«,))2/^'D„(n)-^/^'d"w (D.ll) 







This integral can be calculated only for this special relations between the parameters and is equal: 

MP) = cosh ( . - + ^ r(i + 5 + ¥)r(i + 5-¥) 



X 



r(i - g) 7re-(-/^'+^^/'')r(i + ^^) 
lir(i-^)r(i + | + .f)r(i + f^-f^) ^^'"'^ 



All other normalization integrals for conformal blocks considered in the appendix [E] can be expressed 
in terms of the integrals /n(4n/62, 2P/6, -2/6^) and /„((4n - 2)/62, 2P/6, 
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Appendix E. Integrals for conformal blocks 



Here we list elliptic conformal blocks which have explicit rather simple integral representation. As 
was proposed in section [3] elliptic conformal blocks which have explicit integral representation can be 
labeled by four integer numbers and for k = 1, 2, 3, 4 (where Sk = + ^) 

(E.l) 



(mi,ni) (m4,n4) 



def-,/2 4 b 2 4 b 

■'J^ ' Q _ (2mi+l)fc _ ni ± Q _ (2m4+l)fe _ 714 
2 4 6+26 2 4 b 



where the conformal block in the r.h.s. of (lE.ip is defined as in (ll.lSp with the identification of the 
points X, 0, 1 and oo as follows 

(ml,nl)^^x, (m2,n2)i-^0, (m3,n3)^^oo, (m4,n4)i-^l. (E.2) 

Blocks Hp{. . . \q) are more complicated and we will not consider them in this appendix (see some dis- 
cussions at the end of the section [3]). We give expressions for the conformal blocks up to normalization 
which can be always expressed in terms of F-functions using integrals given in the appendix [Dl To 
reduce the integrals considered in this appendix to normalization integrals we should take into account 
that in the limit q 

Ei{u) sm{u), E2{u) cos(m), Esiu) 1, E^{u) 1, (E.3) 

where the functions Ek{u) are given by flA.131) . Below we give expressions for the conformal blocks 
(lE.ip mostly in cases when only several of numbers m^'s or ra^'s are non-zero. We adopt a notation 
where we will write only the parameters rrikS or n^'s which are non-zero (see below). 

Firstly we consider the situation when only some numbers m^'s are non-zero. The number of 
integrations coinsides in this case with the number of gaps fl3.14p in the potential (12.130 . In the case 
ruk = m and other mj = the conformal block is given by m-dimensional integral 

7r TV 

HUm, = m) ~ Uf[e''''''^E,{uX''^^ J] lEiM^' d^u. (E.4) 

Q Q a=l ^ ' a<a' 

where d"^u = dui...dum- Now we consider the situation when two numbers in (IE. II) are non-zero. 
For the case mi = m, = m and other rrij = we have 

„ » m 

H+{m, = m, = m)^ \[e''''^-{E,Ek{ua)f"^'' \[ \E,Ek{uaa')\-^'' d^u, (E.5) 

a=l a<a' 







where we denote for shortness EiEk{u) == Ei{u)Ek{u). For mi = 0, rrii = m, rrij = m and other 
irik = the conformal block is given by 

i7+(m, = m, = m) ~ f[e''^^^{EiE,{u^)r-'' ft^^'"''] J] \E,E,{u^,,)\-''^' d"^u, (E.6) 

Q Q a=l a<a' 
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Now we consider the less trivial case of non-coinciding integers mi = m, = m — 1 and rrii = rrij = 
with i, j ^ 1 ^ k. The integral in this case will be m-dimensional 

Ht{m, = m,m, = m - 1) ~ j ...j \{e''^'^'^{E,E,{ua)f'-'' { eMu] T H lEiE^Mr'"' d^u. 

Q Q a=l ^ a<a' 

(E.7) 

For the case when all numbers = m conformal block can be derived from eq (lE.4p with /c = 1 by 
substitution 6 ^ 2b. 

Now we consider situation when some numbers n^s in flE.ip are non-zero. In this case the dimension 
of the integral is always equal to = Ylj ^j- the case Uk = n and other rij = the conformal block 
is given by a n-dimensional integral 



In the case ni = n, = n one has 
Htirii =nk =n) 



j...j\{e''-'''--{E^Ek{ua)p n \EiEuM\-^ <f-u. (E.9) 

Q Q a=l a<a' 

We note that the conformal block (]E.9P is dual to the conformal block (]E.5|) (with m substituted by 
2?T, and 6 ^ in (1E.5I) ). In the case Ui = n, rij = n for i,jy^l one has 

1" 2n 1 

HUn^=n,=n) ~ H e^^"^-(i5;ii5;.K))^ (ff^) H l^i^^KaOP^ (E.IO) 

a=l ^ a<a' 







This conformal block is dual to the conformal block (1E.7P (with m substituted by 2n and b ^ b ^). 
We give here examples of non-coinciding numbers rij. For the case rii = n — 1, rik = n with other two 



numbers Uj = one has 



2n-l 

2 



/f+(ni = n-l,nfc = n)~ /.../ J] e''"''""(i^^ii5;fc(Mj)^ J] |EiEfc(M,,Or^ (E.ll) 







a=l a<a' 



This block is dual to flE.5l) with m — *• 2n — 1 and 6 — >■ 6 ^. It is interesting also to consider the case 
rii = n — 1 and n2 = = = n. The integral will be {An — l)-dimensional 



4n-l 



H+{m = 71-1,712 = ns = n, = n)^ ■■■ J] ^''"''""^i H '^i^^-')'"^ (E'^^) 







a=l a<a' 



which is dual to the conformal block (IE. 40 with m = 4n — 1 and A; = 1. In the case ni = n2 = na = 
n4 = n a conformal block is given by the integral which is again dual to the conformal block (1E.4P 
with m = 4n and k = 1. We plan to discuss this intriguing duality elsewhere. 
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We considered conformal blocks with only numbers m^'s either n^'s not equal to nonzero. In the 
mixed case conformal blocks are given by more complicated integrals. For example 



H^[mk = m,nk = n) YYEi[Ua) ^ ) 11 -^iKO \ Ei{u' ,) ) ' 

\{\E,{vac)\-'' n \EiK,)\~^\{{E,{va~v'^,))-^<rv<rv'. (E.13) 







a<c a' <c' 



Conformal block (1E.13P is originated from the integral (1C.12P and hence the same regularization by 



taking the principal value should be applied. 

In this appendix we considered only several examples of conformal blocks which have an integral 
representation. At present integral representations for conformal blocks corresponding to general values 
of integers and in (13.31) are unknown. At the end of this appendix we note that for all cases 
considered here the corresponding correlation functions can be derived in the same way as it was done 
above (see eq (I3.26P ). Namely we should change all elliptic theta functions Ei{u) in the corresponding 
integrands by the functions 

Ekiu, u) = Ek{u)Ek{u)e~^^ (E.14) 

multiply by (Imr)^^, put P = and perform the integrals in the variables u,u over a torus with 
periods vr and ttt. 
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